RESCEU-5/12 



(N 

o 

(N 



Cosmology based on f(R) Gravity admits 1 eV Sterile Neutrinos 
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It is shown that the tension between recent neutrino oscillation experiments, favoring sterile 
neutrinos with mass of order of 1 eV, and cosmological data which impose stringent constraints 
on neutrino masses from the free streaming suppression of density fluctuations, can be resolved in 
models of the present accelerated expansion of the Universe based on f(R) gravity. 
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The possibility of the existence of light sterile neutri- 
nos has been suggested by neutrino oscillation experi- 
ments such as the Liquid Scintillator Neutrino Detector 
(LSDN) y] and the Booster Neutrino Experiment (Mini- 
BooNE) 0]. Recent nuclear reactor experiments also fa- 
vor the additional neutrinos with mass of this range |3| . 
Further analysis shows increasing experimental evidence 
that there may exist one or two species of sterile neutri- 
nos with mass of order of 1 eV [J-Q • 

From the cosmological point of view, recent reanalysis 
of the primordial helium abundance produced at the big 
bang nucleosynthesis (BBN) also favors the existence of 
extra components of radiation [7| . In terms of the effec- 
tive number of neutrinos, which is defined by the total 
energy density of the radiation as 
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with p~ 1 being the energy density of photons, they find 



N, 



off 



3.68 



+0.80 



(2a) or N cS = 3.80 
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-0.70 



(2a) for the 



neutron lifetime r n = 885.4 ± 0.9s or 878.5 ± 0.8s, re- 
spectively. Note that the standard three flavor neu- 
trino species gives N e g = 3.046 Q. Furthermore, the 
cosmic microwave background (CMB) anisotropy obser- 
vations at the small angular scales yield similar value, 
N cS = 4 - 5 0-QJ1. Thus both BBN and CMB suggest 
there exists extra relativistic species, which may be ster- 
ile neutrinos that are expected to be thermalized in the 
early Universe due to mixing [13| . 

If we further incorporate the large-scale-structure 
(LSS) data in the cosmological analysis, however, it turns 
out that in the standard flat A cold-dark-matter (ACDM) 
model the sterile neutrino mass is constrained to be ap- 
preciably smaller than 1 eV fl4T - [l6J to avoid suppression 
of small-scale fluctuations due to free streaming. 

Hence apparently there exists some tension between 
experimental data of neutrino oscillations and cosmolog- 
ical data. But of course they should not be treated on 
equal footing because the latter requires a number of as- 
sumptions on the cosmic evolution. Indeed there have 



been some attempts to make cosmology compatible with 
these experimental data [13, [lj| by adopting a wCDM 
model to treat the equation-of-state parameter of the 
dark energy, w, as an additional fitting parameter, or 
by introducing extra radiation components besides the 
sterile neutrinos. However, the former results in w < —1 
and larger CDM abundance with the cosmic age signifi- 
cantly smaller than the standard value, while the latter 
solution may be in conflict with the aforementioned con- 
straint from BBN. So neither of them is an attractive 
solution. 

In this Letter, we show that extra growth of small- 
scale fluctuations at recent redshifts which occurs in vi- 
able cosmological models of the present accelerated ex- 
pansion of the Universe (in other terms, in models of 
present dark energy) based on f(R) gravity |19rl2l| can 
make eV-mass sterile neutrinos compatible with cosmo- 
logical observations under the proper choice of the func- 
tion /. f(R) gravity is a simple generalization of General 
Relativity (GR) obtained by introducing a phenomeno- 
logical function of the Ricci curvature R, see e.g., the 
recent review [22j |. It represents a special case of more 
general scalar-tensor gravity with the Brans-Dicke pa- 
rameter lobd = 0, and it has an extra scalar degree of 
freedom (or, scalar particle dubbed scalaron). However, 
in contrast to Brans-Dicke gravity, scalaron is massive 
and its rest mass M s depends on R , i.e., on background 
matter density in the regime of small deviations from 
GR. Such models explain the present cosmic acceleration 
without introducing of a cosmological constant, mathe- 
matically this means that /(0) = 0. For the f(R) models 
of present dark energy constructed in [19h21| which sat- 
isfy all existing observational data, the deviation of the 
background evolution from the standard ACDM model 
is small, less than few percent (see e.g., [23[ ), though not 
exactly zero. Their most dramatic difference from the 
standard model appears in the enhancement of the evo- 
lution of matter density perturbations on scales smaller 
than the Compton wavelength of the scalaron field that 
occurs at redshifts of the order of a few depending of the 
scale. This extra growth can compensate the suppression 



due to the free streaming of massive neutrinos and thus 
the upper bound for neutrino mass is relaxed in f(R) 
gravity [24|. We show that the same mechanism works 
in the case of sterile neutrinos, too, to make cosmology 
with them compatible with neutrino experiments. 
f(R) gravity is defined by the action, 
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where S m is the action of the matter content and it is 
assumed to be minimally coupled to gravity (we put H = 
c = 1). If we set f(R) = R — 2A, it reproduces GR with 
a cosmological constant. Instead, for definiteness we use 
the following form [2l| 



f(R) =R + XR 8 
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where n, A, and R s are model parameters. The two of 
them are free parameters and the other one is determined 
by other two and observational data. If we take n and 
A as free parameters, R s is approximately proportional 
to A -1 [23(. Note that n should be taken sufficiently 
large, n > 2, if we want to obtain a noticeable effect for 
the density perturbation enhancement (see below). The 
model (O can describe the accelerated expansion of the 
present Universe and it quickly approaches the ACDM 
model for redshift z > 1 if we take large n and A. Strictly 
speaking, a term proportional to R 2 should be added to 
© to avoid the scalaron mass M s exceeding the Planck 
mass for high, but not too high matter densities in the 
early Universe [2l| , as well as to exclude possible forma- 
tion of an extra weak curvature singularity in the recent 
past which was found in [23, [26[ (still there remains an 
open question what would occur instead of this singular- 
ity j27|)- However, the coefficient of this term (usually 
written as (6M 2 ) -1 where M coincides with the scalaron 
mass M s in the regime when this term dominates other 
non-GR terms in ©) should be very small in order not 
to destroy the standard evolution of the early Universe. 
Namely, either M should be larger than the Hubble pa- 
rameter H at the end of inflation, or this term can drive 
inflation by itself if M « 3 X 10 13 GeV [Hill. With the 
R 2 term, the functional form ([3]) which is responsible for 
the late time acceleration should also be reconsidered and 
changed for R < R s (including the region R < 0) that is 
needed for nonsingular evolution of the model just after 
the end of inflation. We study this point deeply in [301 ■ 
But if we are interested in the dynamics of the late time 
Universe only, we can neglect both the R 2 term and the 
accompanying modification for R < R s . Hence, we use 
Eq. ([3]) as the low energy effective theory for the follow- 
ing. 

The model ^ describes similar background expan- 
sion history with that of the ACDM model. Although 
the equation of state parameter for dark energy makes 
a phantom crossing at z ~ 3 [lj| [23|, [31(, it does not 
change CMB spectrum significantly. On the other hand, 



fluctuations evolve differently. We define the metric per- 
turbation by the following notation, 

ds 2 = -(1 + 2<$>)dt 2 + a 2 (t)(l - 2V)8 ij dx i dx j . (4) 

We can derive the effective gravitational constant and 
the gravitational slip in f(R) gravity by using subhorizon 
limit and quasi-static approximation: 
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Here, A is the gauge-invariant comoving matter pertur- 
bation, and the prime denotes the derivative with respect 
to the Ricci curvature. Thus, in the quasi-GR regime 
when f«l, the effective gravitational constant can be- 
come up to 33% larger, independently of a detailed func- 
tional form of f(R). This is the cause of the enhancement 
of perturbation growth involved. 

As a result of the time and the scale dependences of 
these parameters, evolution of matter density fluctuation 
is different from that in ACDM model, namely, it is en- 
hanced. It promotes formation of LSS [l9l 1 2ll l23l. l32l| . On 
the contrary, the light neutrinos suppress structure for- 
mation by free streaming. Therefore, f(R) modification 
and neutrino masses play opposite roles on the growth 
of perturbation and thus the allowed range for the total 
neutrino mass is relaxed in f(R) gravity, compared with 



24j . We can apply this mechanism for 



the ACDM model 

the case of sterile neutrinos. 

We have carried out Markov Chain Monte Carlo 
(MCMC) analysis for the ACDM model and f(R) grav- 
ity with one or two sterile neutrinos. Practically, we ne- 
glect rest masses of three standard neutrino types (as- 
suming that Yli=i m vi < 0.1 eV) compared to those 
of one or two sterile neutrino types. We have modi- 
fied the MGCAMB (HHl], which provides the evolution 
of the modified growth of matter fluctuation by setting 
G e s(t, k) and r](t, k) as special parameterization, so that 
it allows to implement f(R) gravity by adopting (|7|) and 
©. We do not change the background evolution equa- 
tions, i.e., we keep that in the ACDM model, because 
the difference between the background evolution in the 
viable f(R) model and the ACDM model is not signifi- 
cant (though it is not exactly zero). We used the above 
modified MGCAMB and CosmoMC [H, [Hi to constrain 
the model parameters. The free parameters are the den- 
sity parameter for the dark matter SIdm^ 2 , sound horizon 
angle 8* = \0§r s (z*) / D a{z*) , massive neutrino ratio, /„, 
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FIG. 1: 1 and 2<r contours of the sterile neutrino mass (left) 
and as (right) for the cases with three massless and one mas- 
sive neutrinos in the ACDM model (dashed black) and f(R) 
gravity (solid blue), xls = 3774.1 and 3767.0, respectively. 



FIG. 3: 1 and 2a contours for the sterile neutrino mass and 
f(R) parameter A for the case with three massless and one 
massive neutrinos (left) and three massless and two massive 
neutrinos (right). A > 0.95 to guarantee the stability of the 
future de Sitter stage. 
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FIG. 2: 1 and 2a contours for the sterile neutrino mass (left) 
and as (right) for the cases with three massless and two mas- 
sive neutrinos in the ACDM model (dashed black) and f(R) 
gravity (solid blue), xls = 3788.7 and 3779.3, respectively. 



which is the fraction of the energy density of dark matter 
in the form of massive neutrinos: 



fu = 
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and the amplitude of the f(R) modification, A. Here, 
f^DM means the sum of the contribution from cold dark 
matter and massive neutrinos. The value n — 2 is cho- 
sen because it is the minimal integer value for which the 
scalaron mass M s given by M 2 = 1/3/" (R) in the quasi- 
GR regime is, on one hand, not much higher than the 
Hubble constant Ho if estimated at the present cosmic 
background matter density p m o = 3Q m HQ / (8nG) (see 
the value of the constant Bq below which characterizes it 
quantitatively), and on the other hand, it is already suffi- 
ciently large for matter densities inside the Solar system 
(not speaking about those in laboratory experiments) 
to make scalaron heavy and unobservable even outside 
gravitating bodies similar to dilaton in the string theory. 
Indeed, for the functional form (O, M s ex p™ +1 in the 
quasi-GR regime for p m 3> p m o- Thus, here there is no 
necessity to consider the more subtle chameleon effect 
which can make scalaron unobservable even if it is light 
outside bodies (though heavy inside them). 

Further, in order to have the future stable dc Sitter 



stage, A should be larger than 0.95. For n = 2 and A = 
0.95, the deviation index B = (/"//') (dR/d In H)\ t=tQ 
is not too small nor too large, namely, Bo — 0.21 [23j j . 
This value is in agreement with the upper limit Bq < 0.4 
recently obtained in [37| ■ Contrary, the much more strin- 
gent upper limit obtained from cluster abundance in 
[38( is not applicable for our model because it was ob- 
tained while for a similar functional form of f(R) intro- 
duced in J19l | , but for its parameter value characterizing 
the large- R behaviour which corresponds to n = 0.5 in 
Eq. ©. We fix the other cosmological parameters by 
the mean value constrained by WMAP7 data only. To 
constrain the free parameters, we used the observational 
data of CMB b y W MAP7 @ and matter power spectrum 
by SDSS DR7 [3|. 

Figures [T] and [2] depict the contour plot for the case of 
one and two sterile neutrinos, respectively. In the ACDM 
model, the total neutrino mass is constrained in the sub- 
eV range at 2a. On the other hand, it is allowed to 
take up to the order of 1 eV in f(R) gravity. We also 
find a larger value of as with the mean value about 0.8 
that expresses the enhancement of perturbations as we 
mentioned above. Comparing Figs. [T] and [51 it is seen 
that in the two sterile neutrino case the value of ag is 
slightly smaller, about 5 %. This presents the possibility 
to distinguish these two scenarios when more exact data 
on as will be obtained, e.g., from cluster abundance. 

Figures [3] also suggest f(R) gravity favors the total 
sterile neutrino mass around 1 eV with mild values of 
A with which small-scale enhancement of fluctuations is 
appreciable. 

Moreover, we find that in the presence of the massive 
sterile neutrino f(R) gravity fits the cosmological data 
much better than the ACDM model. Indeed, for the case 
with a single sterile neutrino species the best-fit values 
are m v = 0.860 eV, A = 5.72, and xls = 3767.0 for 
f(R) gravity, and m v = 0.109 eV and xts = 3774 - X for 
the ACDM model. According to the Akaike Information 
Criteria (AIC) [40], if xts improves by 2 or more with 
a new additional fitting parameter, its incorporation is 



justified. In this context, the performance of f(R) gravity 
is very well improving XeS by more than 7. In the case 
with two sterile neutrino species, the improvement of xts 
is even larger, 9.4. 

Note that these are comparisons between the best-fit 
values for which the ACDM model gives a significantly 
smaller values of the sterile neutrino mass. Needless to 
say, the existence and the mass of the sterile neutrino is to 
be determined by ground-based experiments rather than 
cosmology which also depends on other factors that fix 
the evolution of the background and perturbed Universe 
in a complicated manner. 

If, for example, future experiments fix the sterile neu- 
trino mass at, say, m u = 1 cV, then f(R) gravity yields 
xl s = 3767.3 with the best-fit A = 5.97, while the 
ACDM model yields XeS = 3778.0 significantly worse 
than the former. Moreover, ACDM model yields erg = 
0.622 ± 0.005, which is in disagreement with observa- 
tion [ill (39|. f(R) gravity, on the contrary, also works 
well in this respect with erg = 0.76 ± 0.02 for this fixed 



mass of m v . 

Thus we conclude that if a ~ 1 cV sterile neutrino is 
indeed established by ground-based experiments, that is 
in fact favored by a number of experiments now, then cos- 
mological data strongly favor cosmology based on f(R) 
gravity compared to the canonical ACDM model. 
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